Evidence for muon neutrino oscillation in an accelerator-based experiment 
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(Dated: February 7, 2008) 

We present results for z/ M oscillation in the KEK to Kamioka (K2K) long-baseline neutrino os- 
cillation experiment. K2K uses an accelerator-produced beam with a mean energy of 1.3 GeV 
directed at the Super-Kamiokande detector. We observed the energy dependent disappearance of 
which we presume have oscillated to v T . The probability that we would observe these results if 
there is no neutrino oscillation is 0.0050% (4.0a). 

PACS numbers: 14.60.Pq,13.15.+g,25.30.Pt,95.55.Vj 



Recent atmospheric 0, 0, 0, reactor 0, and solar 
neutrino experiments show that the existence of 

neutrino oscillation and non-zero neutrino mass are very 
likely. Measurements of atmospheric neutrino suggest v„ 
to v T oscillation with a mass squared difference (Am 2 ) 
around 2.5 x 10 _3 eV 2 and a mixing angle parameter 
(sin 2 29) that is almost unity 0,0. 

The KEK to Kamioka long-baseline neutrino oscilla- 
tion experiment (K2K) 8, 9] is the first accelerator based 
project to explore neutrino oscillation in the same Am 2 
region as atmospheric neutrinos. The neutrino beam is 
98% Vp, whose direction is monitored every beam spill 
by measuring the profile of muons from the pion de- 
cays. The neutrino beam energy spectrum and profile 
are measured by the near neutrino detectors located 300 
m from the production target. They consists of two de- 
tector sets: a 1 kiloton water Cherenkov detector (1KT) 
and a fine grained detector system. The far detector is 
Super-Kamiokande (SK), a 50 kiloton water Cherenkov 
detector, located 250 km from KEK. 

In this letter, we present evidence for the energy- 
dependent disappearance of v^, which are presumed to 
have oscillated to v T . We observe a distortion of the neu- 
trino energy (E v ) spectrum and a deficit in the total num- 
ber of events. The expectation for these are derived from 
measurements at the near detectors and transformed us- 
ing the energy-dependent ratio of the flux at the far 
and near detectors (F/N ratio). This ratio accounts for 
the difference between the small portion of the beam near 
the center seen by SK and the large section of the beam 
seen by the near detectors. This is calculated using the 
neutrino beam Monte Carlo (MC) simulation and con- 
firmed by measurements of pions from the target . 

Data sample. — We have analyzed data taken from 
June 1999 to February 2004, which corresponds to 8.9 x 
10 19 protons on target (POT). From 1999 to 2001 (called 
K2K-I and SK-I), the inner detector surface of SK had 
11,146 20-inch photo-multiplier tubes (PMTs) covering 
40% of the total area 10]. The fine grained detector 
system was co mpri sed of a scintillating fiber and water 
detector (SciFi) [llj, a lead glass calorimeter, and a muon 
range detector (MRD) 12]. Starting from January 2003 
(K2K-II and SK-II), 19% of the SK inner detector is cov- 
ered using 5182 PMTs, each enclosed in a fiber reinforced 
plastic shell with an acrylic cover. The transparency and 
reflection of these covers in water are 97% and 1% re- 
spectively. The near detector data from this period in- 



clude 2.3 x 10 19 POT without the lead glass (K2K-IIa), 
and then 1.9 x 10 19 POT (K2K-IIb) with a fully-active 
scintillator detector (SciBar) 01 m its place. Adding 
the K2K-II data doubles the statistics compared to the 
previous analysis 0. The neutrino beam direction is 
monitored using neutrino events in the MRD. It is sta- 
ble, within 1 mrad throughout the entire experimental 
period. Also, these events confirm that the energy spec- 
trum is stable. 

The 1KT data alone is used to estimate the expected 
total number of events at SK because the 1KT uses the 
same water target and the uncertainties in the neutrino 
cross section cancel. The event selection and the 25 ton 
fiducial volume are the same as in I9|. We select the 
subset of events in which all the energy is deposited in 
the inner detector (fully contained) and only one, muon- 
like Cherenkov ring is reconstructed (1-ring /it- like events) 
to estimate the E v spectrum along with data from the 
other near detectors. For these events, we reconstruct 
E v by using the measured muon momentum (p^) and 
direction (9^). For the energy spectrum measurement, 
the largest contribution to the systematic uncertainty is 
+2/-3% in the overall energy scale. 

The SciFi detector is made of layers of scintillating 
fibers between aluminum tanks filled with water. The 
fiducial mass is 5.6 tons. The K2K-I analysis includes 
events which reach the MRD and also events in which 
the muon track stops in the lead glass, with momentum 
as low as 400 McV/c, significantly lowering the energy 
threshold compared to [9j . The muon momentum thresh- 
old for K2K-IIa is 550 MeV/c because in this case we 
restrict our analysis to events which have hit at least two 
layers in the MRD to improve the purity of muons. 

The SciBar detector consists of 14,848 extruded scin- 
tillator strips read out by wavelength shifting fibers and 
multi-anode PMTs. Strips with dimensions of 1.3 x 2.5 x 
300 cm 3 are arranged in 64 layers. Each layer consists of 
two planes to measure horizontal and vertical position. 
The scintillator also acts as the neutrino interaction tar- 
get; it is a fully active detector and has high efficiency 
for low momentum particles. Although the target is not 
water, possible differences due to nuclear effects are in- 
cluded in the systematic uncertainty. 

In SciBar, tracks which traverse at least three layers 
(~8 cm) are reconstructed. The reconstruction efficiency 
for an isolated track longer than 10 cm is 99%. In the 
present analysis, we select charged current (CC) events 
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TABLE I: The reconstruction efficiency [%] and purity (in 
parentheses, [%]) for the quasi-eiastic interaction in each sub 
sampie estimated by MC simuiation. 





1-track or 


2-track 
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QE 


non-QE 




1KT 


53 (59) 
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39 (50) 


5 (53) 


2(11) 


46 
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36 (57) 


5 (58) 


2 (12) 


42 


SciBar 


51 (57) 


15 (72) 


4(17) 


70 


SK 


86 (58) 
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by requiring at least one of the tracks start from the 
9.38 ton fiducial volume and extend to the MRD. With 
this requirement, the p M threshold is 450 MeV/c. The 
Pfj, scale uncertainty, p^ resolution, and 9^ resolution are 
2.7%, 80 MeV/c, and 1.6°, respectively. The efficiency 
for a second, short track is lower than that for a muon 
track mainly due to the overlap with the primary track. 
This efficiency smoothly increases from the threshold (8 
cm, corresponding to 450 MeV/c proton) and reaches 
90% at 30 cm (670 MeV/c for proton). 

For SciFi and SciBar, we select events in which one 
or two tracks are reconstructed. For two-track events, 
we use kinematic information to discriminate between 
quasi-elastic (QE) and non-QE interactions. The direc- 
tion of the recoil proton can be predicted from p^ and 9^ 
assuming a QE interaction. If the difference between the 
observed and the predicted direction of the second track 
is within 25°, the event is in the QE enriched sample. 
Events for which this difference is more than 30° (25°) 
for SciFi (SciBar) are put into the non-QE sample. The 
QE efficiency and purity of the samples are estimated 
from the MC simulation and are summarized in Tab. [TJ 

Near detector E v spectrum. - We measure the 
E v spectrum at the near detectors by fitting the two- 
dimensional distributions of p^ versus 9^ with a baseline 
MC expectation . We simultaneously obtain the cross 
section ratio of non-QE to QE interactions (R nqc ) relative 
to our MC simulation. However, we observe a significant 
deficit of forward going muons in all near detector data 
compared to the MC. To avoid a bias due to this, we 
perform the E v fit using only data with 9^ > 20(10) de- 
grees for 1KT (SciFi and SciBar). The x 2 value at the 
best fit is 538.5 for 479 degrees of freedom (DOF). The 
resulting E v spectrum and its error are summarized in 
Tab. ITU while the best fit value of i? nqc is 0.95. 

Muons in the forward direction also correspond to 
events with a low value for the square of the momentum 
transfer (q 2 ), the relevant parameter in the neutrino in- 
teraction models. From inspection of all subsamples, the 
amount of resonant pion production and coherent pion 
production at low q 2 in the MC simulation are possible 
sources of the forward muon deficit. In our MC, we use 
the model for resonant pion by Rein and Sehgal |l4| with 




FIG. 1: The q^ oc distributions for 2-track non-QE samples 
of SciFi (left) and SciBar (right). Open circles, solid lines, 
and hatched areas show data, MC predictions, and CC-QE 
component estimated from MC simulation, respectively. 



axial vector mass of 1.1 GeV/c 2 . For coherent pion, we 
use the model by Rein and Sehgal 115 1 with the cross 
section calculated by Marteau et al. |l6j. Figure shows 
the q 2 distributions calculated from p M and 9^ assuming 
CC-QE kinematics (<z 2 oc ). We modify the MC simulation 
used in the near and the far detector analysis to account 
for the effect of the observed deficit. For resonant pions, 
we suppress the cross section by q 2 /A for q 2 < A and 
leave it unchanged for q 2 > A. From a fit to the SciBar 
2-track non-QE sample, A is 0.10 ± 0.03 (GeV/c) 2 . Al- 
ternatively, if we assume that the source of the low q 2 
deficit is coherent pion production, we find the observed 
distribution is reproduced best with zero coherent pion. 

Considering both possibilities mentioned above, we fit 
the parameter i? nqo again and check the agreement with 
the data. The E u spectrum is kept fixed at the values 
already obtained in the first step, but now we use data at 
all angles. The best fit value for R nqc is 1.02 (1.06) with 
X 2 /DOF of 638.1/609 (667.1/606) when we suppress res- 
onant pion (eliminate the coherent pion). The p^ and 
9^ distributions from all detectors are well reproduced 
for both cases with reasonable x 2 , as shown in Fig. [21 
If we repeat the fit with the E v spectrum free, the re- 
sults are still consistent with the first step. Examining 
these results carefully, we conclude that we cannot iden- 
tify which is the source of the observed deficit in the low 
q 2 region. Because the value of i? n qo changes depending 
on the choice of model, an additional systematic error of 
0.1 is assigned to i? n qo- For the oscillation analysis pre- 
sented in this letter, we choose to suppress the resonance 
production mode in the MC simulation and when we de- 
termine the central value of i? nqc . However, we find that 
the final oscillation results and allowed regions do not 
change if we instead choose to eliminate coherent pion, 
or use our MC without any corrections. 

Oscillation analysis. - Events in SK from the ac- 
celerator are selected based on timing information from 
the global positioning system. The background coming 
from atmospheric neutrinos is estimated to be 2 x 10~ 3 
events. For K2K-I+II there are 107 events in the 22.5 
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TABLE II: The E v spectrum fit results. $nd is the best 
fit value of relative flux for each E v bin to the 1.0-1.5 GeV 
bin. The percentages of uncertainties in <I?nd, F/N ratio, and 
reconstruction efficiencies of SK-I and SK-II are also shown. 
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FIG. 2: A selection of muon momentum (p M ) and direc- 
tion (# M ) distributions: (a) the distribution of 1KT fully 
contained 1-ring /i-like sample, (b) 1KT 9^ for the same sam- 
ple, (c) SciFi for 2-track QE sample, and (d) SciBar # M 
for 2-track nonQE sample. Open circles represent data, while 
histograms are MC predictions using the best fit E v spectrum 
and suppression of the resonant pion production. 



kiloton fiducial volume that are fully contained, have no 
energy seen in the outer detector, and have at least 30 
MeV deposited in the inner detector. The expected num- 
ber of fully contained events at SK without oscillation is 
151±io(syst). The major contributions to the errors come 
from the uncertainties in the far to near ratio (5.1%) and 
the normalization (5.1%); the latter is dominated by the 
uncertainty in the fiducial volumes due to the vertex re- 
construction at both 1KT and SK. 

We reconstruct the neutrino energy (E^ oc ), assuming 
CC-QE kinematics, from and 9^ for the 57 events in 
the 1-ring /i-like subset of the SK data. With these we 
measure the energy spectrum distortion caused by neu- 
trino oscillation. The detector systematics of SK-I and 
SK-II are slightly different because of the change in the 
number of inner detector PMTs. In the oscillation anal- 
ysis based on the energy spectrum, the main contribu- 
tion to the systematic error is the energy scale uncer- 
tainty: 2.0% for SK-I and 2.1% for SK-II. Uncertainties 



for the ring counting and particle identification are esti- 
mated using the atmospheric neutrino data sample and 
MC simulation. The differences between the K2K and 
atmospheric neutrino fluxes are also taken into account. 

A two-flavor neutrino oscillation analysis, with dis- 
appearance, is performed using a maximum-likelihood 
method. The oscillation parameters, (sin 2 29, Am 2 ), are 
estimated by maximizing the product of the likelihood 
for the observed number of FC events (£ num ) and that 
for the shape of the E™ c spectrum (£ s ha P e)- The prob- 
ability density function (PDF) for £ num is the Poisson 
probability for the expected number of events. The PDF 
for £ s h a pe is the expected E™ c distribution at SK, which 
is estimated from the MC simulation. The PDFs are de- 
fined for K2K-I and K2K-II separately. The systematic 
uncertainties due to the following sources are taken into 
account in the PDFs: the E v spectrum measured by the 
near detectors, the far to near ratio, the reconstruction 
efficiency and absolute energy scale of SK, the ratio of 
neutral current to charged current QE cross section, the 
ratio of CC non-QE to CC-QE cross section and the over- 
all normalization. The systematic uncertainties modify 
the expected distributions, and each is assumed to follow 
a Gaussian distribution |7j. A constraint term (£ S yst) is 
multiplied with the likelihood for each of these system- 
atics, and £ num x £shapc x Asyst is maximized during the 
fit. The total number of parameters varied in the fit is 
thirty- four. 




E v reo (GeV) 

FIG. 3: The reconstructed E v distribution for the SK 1-ring 
/i-like sample. Points with error bars are data. The solid line 
is the best fit spectrum. The dashed line is the expected spec- 
trum without oscillation. These histograms are normalized by 
the number of events observed (57). 

The best fit point within the physical region is (sin 2 29, 
Am 2 ) = (1.0,2.8 x 10~ 3 eV 2 ). The expected number of 
events at this point is 103.8, which agrees well with the 
107 observed. The best fit E v distribution is shown with 
the data in Fig. [3] The consistency between the observed 
and fit E v distributions is checked using a Kolmogorov- 
Smirnov (KS) test. For the best fit parameters, the KS 
probability is 36%, while that for the no-oscillation hy- 
pothesis is 0.08%. The highest likelihood is at a point 
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(1.5, 2.2 x 10 -3 eV 2 ) which is outside of the physical re- 
gion. The probability that we would get sin 2 29 > 1.5 if 
the true parameters are our best ht physical parameters 
is 13%, based on MC virtual experiments. For the rest 
of this letter we refer only to the physical region best fit. 
The fit results for all the systematic parameters are rea- 
sonable. The fits for the K2K-I and K2K-II sub-samples 
are consistent with the result for the whole sample. 




0.8 

sin 2 29 

FIG. 4: Allowed regions of oscillation parameters. Dashed, 
solid and dot-dashed lines are 68.4%, 90% and 99% C.L. con- 
tours, respectively. 

The possibility that the observations are due to a sta- 
tistical fluctuation instead of neutrino oscillation is es- 
timated by computing the likelihood ratio of the no- 
oscillation case to the best fit point. If there is no os- 
cillation, the probability of this result is 0.0050%(4.0<t). 
When only normalization (shape) information is used, 
the probability is 0.26% (0.74%). Allowed regions for 
the oscillation parameters are evaluated by calculating 
the likelihood ratio of each point to the best fit point 
and are drawn in Fig. 01 The 90% C.L. contour crosses 
the sin 2 29 = 1 axis at Am 2 = 1.9 and 3.6 x 10~ 3 eV 2 . 
The oscillation parameters from the E u spectrum distor- 
tion alone, or the total event analysis alone also agree. 

In conclusion, using accelerator produced neutrinos, 
we see the same neutrino oscillation discovered with at- 
mospheric neutrino measurements. This result is based 
on data from 1999 to 2004, a total of 8.9 x 10 19 POT. 
The observed number of events and energy spectrum of 
neutrinos at SK are consistent with neutrino oscillation. 
The probability that we would see this result if there was 
no oscillation is 0.0050%(4.0u). The allowed regions of 
the oscillation parameters from the K2K experiment are 



consistent with the atmospheric measurements. 
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